Mutual Impedance of Grounded Wires Lying On or Above 
the Surface of the Earth * 

By RONALD M. FOSTER 

This paper presents a formula for the mutual impedance of any insulated 
wires of negligible diameter lying in horizontal planes above the surface of 
the earth and grounded by vertical wires at their four end-points. The 
formula holds for frequencies which are not too high to allow all displace- 
ment currents to be neglected. Tables and curves are given to facilitate 
numerical computation by means of the formula. 

In the expansion of this formula for low frequencies and for any heights 
the first two terms give the direct-current mutual impedance; the third 
term is independent of the heights, thus being identically the same as that 
previously ffcund for wires on the surface. The mutual impedance for 
wires at any heights H and h, with separations large in comparison with 
these heights, is found to be approximately equal to the mutual impedance 
for wires on the surface multiplied by the complex factor [1 + T(H + ft)], 
where T is the propagation constant in the earth. 

THE formula established in a previous paper l for the mutual 
impedance of any grounded thin wires lying on the surface of 
the earth has now been extended to include wires lying in horizontal 
planes above the surface of the earth and grounded by vertical wires 
at their four end-points. As before, we assume the earth to be flat, 
semi-infinite in extent, of negligible capacitivity, of uniform resistivity 
p, and of inductivity v equal to that of free space. The air is also 
assumed to be of negligible capacitivity and of inductivity v equal to 
that of free space. All displacement currents are thus neglected both 
in the earth and in the air; this is the assumption which is ordinarily 
made as a first approximation at power frequencies for the shorter 
transmission lines. 

By the same general method of derivation as before, the extended 
formula is found to be: 

Za . Jf [ WL*> + cos . M(r,II,h) } dSds, (A) 

where 

P(r,H,h) = Po(f) + Pi(r,H + h) - P 2 (r,\H - h\), 

* A brief report of the principal theoretical result obtained in this paper was 
recently published in the Physical Review (2), 41, 536-537 (August 15, 1932). An 
error in the formula for N {r), as printed there, should be corrected: in the denomi- 
nator of the fraction, r 2 should be r 3 . 

1 R M Foster, "Mutual Impedance of Grounded Wires Lying on the surface ot 
the Earth," Bell System Technical Journal, 10, 408-419 (July, 1931); see also Bulletin 
of the American Mathematical Society, 36, 367-368 (May, 1930). 
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M(r,H,h) = M (r) + Mi(r,H + h) - M 2 (r,\H - h\), 
Po(r) =/-- 
„, , iar r° I s 1 - e->[(S+ r ! )>« -Mil ,, sj 

P,M) = *£ [<* log (f ' ± ^'" ± < - fr. + *)» + r ] , 

«*-*-fcfg-«->[ Siff;: ]^fc • 

The integrations in the iterated integral are extended over the two 
wires 5 and 5, lying in planes at heights H and h, respectively. The 
elements dS and ds are separated by the horizontal distance r and 
include the angle e between their directions. The propagation con- 
stant of plane electromagnetic waves in the earth, varying with the 
time as e iul , is T, which equals (icov/p) 112 . All distances are measured 
in meters, Z ]2 in ohms and p in meter-ohms; v has the value 1.256 
X 10 -6 henries per meter; w is equal to 2tt times the frequency; J is 
the Bessel function of order zero. The derivation of the formula is 
outlined in the latter part of this paper. 

The functions P and M are divided into three parts: first, P and 
M , which are functions only of the horizontal distance r; secondly, 
Pi and Mi, which are functions of r and of the sum of the two heights 
77 and h; and thirdly, P 2 and M 2 , which are functions of r and of the 
numerical difference of the two heights. These three parts are 
arranged in the order of relative importance when the heights are 
reasonably small. For zero heights, the functions P and M reduce to 
P and il/ , which are the values previously obtained for wires on the 
surface. For small values of the heights, Pi and Mi are of the order 
of magnitude of the sum of the heights, whereas P 2 and M 2 are of the 
order of magnitude of the square of the difference. 

For some purposes it is convenient to transform formula (A) into 
the alternative expression : 

Zl2 = ff [ * P dSds h) + C ° S € M{rJI ' h) ] dSds > ^ 
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where 
P(r,H,h) = P (r) + P°(r,H,h) + P 3 (r,II + h), 
M(r,H,h) = M°(r,H,h) + M 3 (r,H + h), 

no, rM^ ^f„, rr 2 + (# + ft) 2 ] 1 ' 2 + // + /* 
F>(r,/J,fc) = ^ ( II log [fi + (H _ h)2JI2 + n _ h 



LI . [r 2 + (H + ft) 2 ]" 2 + H + ft 
-I- " log ^ + (// _ ^ )2 -j 1/2 _ H + ft 

+ [r 2 + (/I - ft) 2 ] 1 ' 2 - l> 2 + (H + /0 2 ] 1 ' 2 } 

f ,(M) = ^1" m[(m 2 +"rr 2 + mJ /oMt/M ' 

M°(r,H,h) = ^ | [r 2 + {H _ A) ij/i - O 2 + (// + /0 2 ] 1 ' 2 } ' 



M 3 (r 



*-ar 



M<r 



( M 2 + r 2 ) 1 ' 2 + /* 



Jo(rp)dp. 



The functions P and M are again divided into three parts: first, 
P , the term giving the direct-current mutual resistance; secondly, 
P° and M°, terms giving the mutual impedance on the assumption of 
a perfectly conducting earth; and thirdly, P 3 and M 3 , the correction 
terms for the finite conductivity of the earth. The P° and M° terms 
thus give iuv/Stt times the mutual Neumann integral of the two com- 
plete circuits formed from the actual wire circuits by adding to them 
their reflections in the surface of the earth. 

For small values of I\ the P 3 and M a terms can be expanded as 
follows : 



P 3{r ,s) = ^ { - s log r - s log [(r 2 + * 2 ) 1/2 + *] - r 

+ (r 2 + s 2 ) 1/2 + [2 log 2 + -Kl) + J> 

+ i^r - &s(3r 2 - 2s 2 ) r 2 log r + • • ■ 



^(r,5) = 7- 



1 



[ (r 2 + s 2 ) 1 ' 2 



-§r - j^rMog r + 



(i) 



By means of these expansions, the complete P and M functions, as 
given by formula (5) , can be put into the form : 
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P(r,H,h) = ^ + ~ { - //log {[r* + (H - h)*J'* + H - h\ 



- h log { [r 2 + (// - hyj' 2 -H+h 
+ [r 2 + (// - h)*J* - r 

+ F(H,h,r) + O(r 2 io g r)}, 



(2) 



M(r,H,H) = £ { p, - (/ / _ /;)2]1/2 - | r + 0(r 2 log r) } ■ 

The function F(H,h,T) is of no consequence, since it does not 
involve r; it contributes nothing to the value of the impedance. The 
remaining terms are infinitesimals of order (T 2 log T) for infinitesimal 
values of T; they are thus of higher order than V itself. 

By means of equation (2) we can now show that the first three terms 
in the expansion of Z X2 for low frequencies and for any heights are 
given by 

7 - -E- ( 1 * 1 A_ 1 \ _L i0)V AT 

Zu ~2r\Aa Ab~B^ + Bb) + 17 ^«">«—» 

1/2 

^15a6cos0+ •••, (3) 

where N^s-ic)( S - c ) is the mutual Neumann integral between the two 
circuits formed by the wires 5 and s, lying in planes at heights // and h 
above the earth, grounded by vertical wires at their four end-points, 
and with earth returns, — the four grounding points on the surface of 
the earth being A,B and a,b, respectively. The angle between the 
straight lines AB and ab is designated by 6. N {S -BH—e) is equal to 
Ns$, the mutual Neumann integral between the two wires S and s, 
augmented by terms which depend only on the arithmetical distances 
between eight points, — the four end-points and the four grounding 
points. 

The first two terms in the expansion (3) are precisely the direct- 
current mutual impedance as given ten years ago by G. A. Campbell. 2 
The third term is independent of the heights of the wires; it is thus 
identically the same as the third term previously found for wires on 
the surface. 

The leading term in the expansion of Z 12 for a long straight wire 5 
and any wire 5 located near the midpoint of S, for any heights, is 



/( 



iup . [.v 2 + (// + h)*y* 



2tt 8 [x 2 + (i/-/;) 2 ] 1 ' 2 

i I v 2 i — n?>i/9 i — cos Xfidfi } cos eds, (4) 

2 G. A. Campbell, " Mutual Impedances of Grounded Circuits," Bell System Tech- 
nical Journal, 2, (no. 4), 1-30 (October, 1923). 
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x being the positive horizontal distance from ds to S, and e the angle 
between ds and S. 

This result is derived immediately from formula (B) upon assuming 
5 to be doubly infinite and then integrating over its entire length. 
The first part of the expression comes from the M° function, the second 
part from the M 3 function. The other functions contribute nothing 
to the leading term in the expansion. 

The expression enclosed in braces in (4) is the mutual impedance 
gradient parallel to an infinite wire at a positive horizontal distance x 
from the wire. It agrees with the results published independently by 
F. Pollaczek, 3 J. R. Carson, 4 and G. Haberland. 6 Pollaczek has also 
investigated the case of two grounded circuits of finite length, with 
certain modifications. 6 

For purposes of computation, however, formula (A) is better, in 
general, than formula (JB). A distinct improvement is effected by 
multiplying all distances which occur in (A) by the attenuation con- 
stant, that is, by (o>j>/2p) 1/2 . The numerical value thus obtained for 
any one distance is indicated by a prime accent on the corresponding 
letter. We then find the mutual impedance expressed in the following 
form: 

Zl , m isg^ff [ ^g^ + cos . WM ] «*\ (C) 

where 

Q(t',ir,h') = Q (r') + Qi{r',H' + h') - Q,{r',\II' - h'\), 
N(r',H',h') = iVo(f') + N 1 (r',H' + h') - N 2 (r',\H' - h'\), 

Qo(r') =p> 

3 F. Pollaczek, " Cber das Feld einer unendlich langen wechselstromdurchflossenen 
Einfachleitung," ' Elektrische Nachrichten-technik , 3, 339-359 (September, 1926). 

4 J. R. Carson, "Wave Propagation in Overhead Wires with Ground Return," 
Bell System Technical Journal, 5, 539-554 (October, 1926). 

6 G. Haberland, "Theorie der Leitung von Wechselstrom durch die Erde," 
Zeitschrift fur angewandte Mathematik und Mechanik, 6, 366-379 (October, 1926). 

6 F Pollaczek, "Gegenseitige Induktion zwischen Wechselstromfreileitungen von 
endlicher Lange," Annalen der Physik (4), 87, 965-999 (December, 1928). His 
assumptions regarding conditions at the ground connections seem to depart con- 
siderably from the conditions assumed in the present paper, and moreover his results 
are not expressed in convenient form for direct comparison with the formula given 
above for Z«. 
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tfoO-o =4n - [i + a + *v>-< i+ '>'}, 



/3 

r 



N 



(V + 2t)" 2 - n 
(/x 2 + 2*) 1 ' 2 + „ 



J (r'ti)dn, 



MS) = • ("" (i - e-v) r 

the prime accent applied to any length L indicating the corresponding 
modified length V = (cop/2p) ll2 L. 

As in formula (A), the six constituent functions involved in formula 
(C) are arranged in order of importance: first, Q Q and N , functions 
only of the modified horizontal distance r'; secondly, Qi and iVi, 
functions of r' and of the sum of the two modified heights IV and h'\ 
and thirdly, Q 2 and N 2 , functions of r' and of the numerical difference 
of the two modified heights. 

To assist in the numerical application of this formula, a table of 
values of the real and imaginary parts of N has been computed, for 
all values of r' from to 10, in steps of 0.1. Beyond this range, the 
function is practically equal to the leading term in its asymptotic 
expansion, namely, l/r' 3 . These computed values are also shown 
graphically in Fig. 1. The imaginary part changes sign at approxi- 
mately r' = 3.8, and again at 7.0, oscillating for increasing values of r', 
although approaching zero very rapidly indeed. 

The real and imaginary parts of the functions Qi(r',s') and Ni(r',s') 
are shown in Figs. 2, 3, 4, and 5, for the range of r' from to 10, and 
for the set of values of s' from to 0.2 in steps of 0.02. It is believed 
that this will cover the range of heights likely to be encountered in 
ordinary problems. To cover this range adequately it was necessary 
to show portions of the N x curves with the horizontal scale enlarged 
two and a half times, and with a greatly reduced vertical scale, in 
Figs. 4- A and 5-^4. 

For actual computation, the Q 2 and N 2 functions are already ex- 
pressed in (A) in convenient, closed form, but for purposes of com- 
parison with Qi and N u the corresponding values of Q 2 and N 2 are 
shown in Figs. 6 and 7, which are drawn to the same scales as Figs. 
3 and 5. 

Tables II and III give the corresponding numerical values of Qi and 
Ni for the range of r' from to 1, in steps of 0.1, as well as the values 
for 1.5 and 2. 
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TABLE I 
Real and Imaginary Parts of N (r') 



r' 


Real 


I mag. 


r' 


Real 


Imag. 





0.66667 


00 


5.0 


0.0081667 


-0.00038659 


0.1 


0.61800 


9.33461 


5.1 


0.0076496 


-0.00034816 


0.2 


0.57199 


4.33824 


5.2 


0.0071782 


-0.00031048 


0.3 


0.52860 


2.67724 


5.3 


0.0067477 


-0.00027431 


0.4 


0.48778 


1.85135 


5.4 


0.0063538 


-0.00024017 


0.5 


0.44949 


1.36031 


5.5 


0.0059927 


-0.00020839 


0.6 


0.41363 


1.03722 


5.6 


0.0056609 


-0.00017918 


0.7 


0.38014 


0.81043 


5.7 


0.0053554 


-0.00015262 


0.8 


0.34892 


0.64405 


5.8 


0.0050736 


-0.00012872 


0.9 


0.31987 


0.51804 


5.9 


0.0048131 


-0.00010740 


1.0 


0.29291 


0.42035 


6.0 


0.0045717 


-0.000088557 


1.1 


0.26792 


0.34327 


6.1 


0.0043477 


-O.OOO072O47 


1.2 


0.24480 


0.28161 


6.2 


0.0041392 


-0.000057706 


1.3 


0.22346 


0.23177 


6.3 


0.0039449 


-0.000045358 


1.4 


0.20379 


0.19116 


6.4 


0.0037634 


-0.000034822 


1.5 


0.18568 


0.15785 


6.5 


0.0035936 


-0.000025919 


1.6 


0.16905 


0.13041 


6.6 


0.0034344 


-0.000018472 


1.7 


0.15379 


0.10770 


6.7 


0.0032850 


-0.000012315 


1.8 


0.13981 


0.088878 


6.8 


0.0031444 


-0.0000072892 


1.9 


0.12703 


0.073237 


6.9 


O.OO3012O 


-0.0000032482 


2.0 


0.11535 


0.060227 


7.0 


0.0028872 


-0.0000000570 


2.1 


0.10470 


0.049402 


7.1 


0.0027693 


0.0000024076 


2.2 


0.095002 


0.040395 


7.2 


0.0026578 


0.0000042564 


2.3 


0.086174 


0.032905 


7.3 


0.0025522 


0.0000055886 


2.4 


0.078152 


0.026685 


7.4 


0.0024522 


0.0000064921 


2.5 


0.070871 


0.021526 


7.5 


0.0023573 


0.0000070444 


2.6 


0.064268 


0.017257 


7.6 


0.0022672 


0.0000073128 


2.7 


0.058287 


0.013734 


7.7 


0.0021816 


0.0000073559 


2.8 


0.052874 


0.010835 


7.8 


0.0021001 


0.0000072236 


2.9 


0.047980 


0.0084577 


7.9 


0.0020226 


0.0000069586 


3.0 


0.043558 


0.0065174 


8.0 


0.0019488 


0.0000065967 


3.1 


0.039567 


0.0049415 


8.1 


0.0018785 


0.0000061678 


3.2 


0.035966 


0.0036689 


8.2 


0.0018114 


0.0000056965 


3.3 


0.032719 


0.0026483 


8.3 


0.0017474 


0.0000052028 


3.4 


0.029792 


0.0018364 


8.4 


0.0016863 


0.0000047027 


3.5 


0.027156 


0.0011967 


8.5 


0.0016280 


0.0000042086 


3.6 


0.024782 


0.00069852 


8.6 


0.0015722 


0.0000037302 


3.7 


0.022645 


0.00031616 


8.7 


0.0015190 


0.0000032745 


3.8 


0.020720 


0.00002803 


8.8 


0.0014680 


0.0000028466 


3.9 


0.018987 


-0.00018390 


8.9 


0.0014193 


0.0000024498 


4.0 


0.017427 


-0.00033467 


9.0 


0.0013727 


0.0000020858 


4.1 


0.016021 


-0.00043678 


9.1 


0.0013280 


0.OOO0O17555- 


4.2 


0.014754 


-0.00050056 


9.2 


0.0012852 


0.0000014587 


4.3 


0.013612 


-0.00053454 


9.3 


0.0012443 


0.0000011945 


4"4 


0.012582 


-0.00054572 


9.4 


0.0012050 


0.00000096159 


4.5 


0.011652 


-0.00053979 


9.5 


0.0011673 


0.00000075815 


4.6 


0.010811 


-0.00052138 


9.6 


0.0011312 


0.00000058219 


4.7 


0.010050 


-0.00049420 


9.7 


0.0010966 


0.00000043156 


4.8 


0.0093603 


-0.00046121 


9.8 


0.0010633 


0.00000030402 


4.9 


0.0087349 


-0.00042473 


9.9 


0.0010313 


0.00000019732 


5.0 


0.0081667 


-0.00038659 


10.0 


0.0010007 


0.00000010925 
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Fig. 1 — Real and imaginary parts of N u (r'). 



272 



BELL SYSTEM TECHNICAL JOURNAL 




a 9 10 



Fig. 2— Real part of Qi(r', s'). 
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Fig. 3 — Imaginary part of Qi{r', s'). 
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Fig. 4— Real part of Ni(r', s'). 
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Fig. 4-A — Real part of N[(r', s'), enlarged horizontal scale. 
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Fig. 5 — Imaginary part of Ni(r', s'). 
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Fig. 5-A— Imaginary part of Ni(r', s'), enlarged horizontal scale. 
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Fig. 6 — Imaginary part of Q«(r', d'). 
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Fig. 7 — Imaginary part of N?(r', d'). 
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These tabulated values were computed from the corresponding con- 
vergent series, the first few terms of which are: 

<M»V) = W - h* + • • • 

+ i{- s' log r' + [f log 2 + *(1) + h¥ 



+ h' 2 + 



NxVJ) = is' log [(r' 2 + s' 2 ) 1 ' 2 + s'] 

- ici log 2 + *(d - iy 

+ v - w* + * ,2 ) I/2 - tV' 2 + ••• 



For values of r* greater than 2, sufficient accuracy for ordinary 
purposes is obtained by using the first two terms in the expansions in 
terms of s': 



(5) 



Qi(r',s') = s'QW) + **Oi«(fO + ' 
jV l( y >5 ') = s'Ni^i.r') + s' 2 iV 2) (r') + 



(6) 



where 



Qi«(rO = i[/o(«)^o(«) + Ji(*)£i(«)]. 



Qi™(/) = 



l 



8m 3 



[(1 - 2« 2 ) - (1+ 2«)c- 2 «], 



(7) 



jVifl^r') = ^[1 - 2uIi(u)K (u) - 2h{u)K 1 {u)~\, 

N^\r') = i-t^[(9 - 2m 2 ) - (9 + 18m + 16m 2 + 8w 3 )e- 2 "], 

u = *(l + *>'■ 

For actual computation we note that 

<2i (2) (r') -*[p- WrO]- 

The real and imaginary parts of these four functions are shown in 
Figs. 8, 9, 10, and 11. The dominating terms in the asymptotic 
expansions of Qi and N% are thus given by those of <2x (1) and N^ l K 
For large values of r', <2i (1) approaches zero as (1 + i)/r', and 7Vi (I) 
as(l+*)/r' 8 . 
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For very large values of r' it is convenient to express the functions 
as follows: 



Go(r') + 0i(r 



',*') = Oo(r') [ 



+ <2o(r') + 



iVo(r') + /V^r'.s') = iV (r') [l + ^y^' + •••] 



(8) 



The real and imaginary parts of these ratios of functions — the coeffi- 
cients of s' in the above expansions — are shown in Figs. 12 and 13. 
We note that each of these ratios approaches the value (1 + i) as f' 
increases without limit. Hence, as a rough approximation, we may 
say that the mutual impedance for wires at heights 77 and h, with 
separations large in comparison with these heights, is equal to the 
impedance for wires at zero heights multiplied by the factor: 



1 + (1 + *)(//' + h') = 1 + T(H + *). 



(9) 



The mutual impedance formula (A) was originally derived from 
first principles, following the method used in the previous paper for 



1.2 



1.0 



a. 0.7 



0.6 
0.5 



| 0.4 
< 



\ 
















































































\ 


\ IMAGINARY 
\ PART 














\ 


V \ 




















\ 




















\ 




















REAL^X 

PART N 


V 











































































4 5 6 

VALUES OF P' 



Fig. 8 — Real and imaginary parts of Qi (l) (r'). 
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Fig. 9 — Real and imaginary parts of Qi <2) (r'). 
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Fig. 10 — Real and imaginary parts of iW l) (r'). 
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wires on the surface. A brief outline of this derivation is given here. 
We first find the formulae for the components of the electric field due 
to a current flowing in a straight wire of length 2a parallel to the 
surface of the earth and at the height H above it, assuming the air 
to be replaced by a medium of finite resistivity p\. This part of the 
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Fig. 11 — Real and imaginary parts of Ni m (r'). 

derivation follows closely the work involved in the previous case of 
wires on the surface. We next find the electric field due to a current 
in a vertical wire extending from the surface of the earth up to the 
height II in the assumed medium. This part of the derivation is 
simpler since there is circular symmetry. Upon combining these two 
results, we obtain the field due to a current flowing through three 
sides of a rectangular circuit beginning and ending at the surface of 
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the earth, extending up to the height H, and of width 2a. We can 
now allow p\ to become infinite, corresponding to the assumptions of 
our problem, since this circuit is completed through the earth. Upon 
allowing a to approach zero, such that 2a = dS, we find the field 
corresponding to a rectangle of infinitesimal width. We then take the 
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Q\ a) (r') 
Fig. 12 — Real and imaginary parts of n , ,. • 

vol.'' ) 

integral of this expression around a similar circuit consisting of a 
horizontal element of wire of length ds at the height h, grounded by 
wires at its end-points. Upon making various algebraic simplifica- 
tions, we finally obtain the mutual impedance as given by formula (A). 
It is perhaps more convenient to derive this formula from results 
obtained by H. von Hoerschelmann, 7 again following the method 

7 H. von Hoerschelmann, "Uber die Wirkungsweise des geknickten Marconischen 
Senders der drahtlosen Telegraphic" Jahrbuch der drahtlosen Telegraphic und 
Telephonic 5, 14-34, 188-211 (September, November, 1911). 
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employed in the previous paper in a similar derivation for wires on 
the surface. For our present problem we use his formulae for the 
Hertzian vectors due to horizontal and vertical electric antennae above 
the surface of the earth. It is important, at first, to retain a non- 
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N\ w (r') 
Fig. 13 — Real and imaginary parts of „ ... • 

vanishing value of the capacitivity of the air. From these formulae we 
obtain the vector II, due to a current flowing through a horizontal 
element of wire of length dS at height H above the surface of the 
earth, grounded by vertical wires at its end-points. Next, we obtain 
the electric field E in the air by the relation : 
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E = grad div n - IVn, (10) 

where T t is the propagation constant in the air. We can now allow 
Ti to vanish, thus obtaining the expression for the field corresponding 
to the assumptions of our problem. We then proceed as before to 
find the expression for the mutual impedance. 

I am greatly indebted to my colleagues, Dr. Marion C. Gray and 
Miss Helen M. Kammerer, for much valuable assistance in the 
preparation of this paper, particularly in the compilation of the tables 
and curves. 



